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Abstract 

DNA sequences from a novel nuclear marker (2nd intron of the vertebrate alcohol dehydro¬ 
genase gene, Adh 1-12) and the mitochondrial cytochrome-/? gene ( Cyt-b ) were examined in 41 
species and six genera of rodents (one Holochilus, 13 Neotoma , six Oryzomys , one Ototylomys, 
19 Peromyscus , and one Tylomys ). The Adh 1 -12 dataset was characterized by having a significant 
level of phylogenetic signal (P< 0.01) and a low level of homoplasy (consistency index = 0.78, 
retention index = 0.89). Based on parsimony, maximum likelihood, and Bayesian analyses, the 
Adh\-\2 results were congruent with results from Cyt-b and other types of data, and generally 
were accompanied by high statistical support. All phylogenetic relationships at and above the 
genus level were strongly supported. Furthermore, compared with Cyt-b , the Adh 1-12 dataset 
provided supported resolution to more species-level relationships within Neotoma , Peromyscus , 
and Oryzomys , despite fewer overall characters (489 Adh 1-12; 1,143 Cyt-b ) and fewer parsimony 
informative characters (126 Adh 1-12; 464 Cyt-b). Overall, results of this study indicated that 
Adh\-\2 DNA sequences are useful for addressing phylogenetic relationships within and among 
the Sigmodontinae andNeotominae. The slower rate of molecular evolution observed \nAdh\- 
12 sequences, coupled with low levels of homoplasy and strength of resolution for supraspecific 
relationships, indicated that this marker may be useful for resolving phylogenetic relationships 
at low to intermediate taxonomic levels in mammals. 
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Introduction 


The use of nucleotide sequence data has greatly 
enhanced our knowledge of mammalian systematics. In 
mammals, the majority of molecular-based phylogenies 
have been estimated by analysis of maternally inherited, 
mitochondrial DNA (mtDNA) sequences, especially 
protein coding genes. Advantages in using mtDNA 
sequences in phylogeny reconstruction include: tech¬ 
nical ease in DNA isolation and polymerase chain 
reaction amplification, manageable sequence lengths 
(generally <1,200 base pairs, bp), ease of alignment, 
availability of'“universal” primers, differential rates of 
molecular evolution in specific regions of the molecule, 
significant phylogenetic signal, and substantial data¬ 
bases (GenBank, EMBL, etc.) of comparative material. 
However, there are drawbacks to mtDNA that can result 
in a misrepresentation of the true phylogeny, includ¬ 
ing retention of ancestral polymorphisms, incomplete 
lineage sorting, differential rates of evolution among 
lineages, and hybridization (Avise 1994; Hillis et al. 
1996; Prychitko and Moore 2000). Likewise, base 
composition of mtDNA protein coding genes (excess of 
C and A at the 3rd positions - Anderson et al. 1981; Roe 
etal. 1985; Desjardins and Morais 1990; Prychitko and 
Moore 2000; Ballard et al. 2002) may result in codon 
bias and homoplasy (Prychitko and Moore 2000) lead¬ 
ing to a biased reconstruction of ancestral sequences 
(Collins et al. 1994; Lockhart et al. 1994; Perna and 
Kocher 1995). Although some of these drawbacks 
can be addressed, nuclear DNA sequences are highly 
desirable as an alternative dataset to independently 
test mtDNA hypotheses and to essentially counteract 
criticisms of mtDNA. 

Historically, nuclear DNA sequences have been 
used to address phylogenetic questions pertaining to 
relationships above the genus level (e.g., rRNA genes 
- Gouy and Li 1989; Perasso et al. 1989), principally 
because nuclear sequence markers with levels of varia¬ 
tion suitable for resolving relationships among closely 
related species or genera are scarcely known and have 
been difficult to discover. Intron markers are an obvi¬ 
ous and likely source of nuclear variability suitable 
for low- or intermediate-level phylogenetics. With the 
exception of some conserved areas, such as regulatory 
sites (Jackson and Hoffmann 1994), most nucleotide 
positions within introns are adaptively neutral and 


independently distributed and generally are evolv¬ 
ing faster than those within exons, thus increasing 
their potential in population genetics or phylogenetic 
analysis of closely related species (Ballard et al. 2002; 
Prychitko and Moore 2000; Slade et al. 1994; Wick- 
liffe et al. 2003). In addition, introns of various sizes 
are relatively abundant in the nuclear genome, and 
usually primer design in conserved areas of adjacent 
exons enable efficient amplification across taxonomic 
boundaries. In vertebrates, several introns have been 
shown to be useful in studies at or below the genus 
level (Prychitko and Moore 1997, 2000; Oakley and 
Phillips 1999; Lavoue et al. 2003; Reeder and Bradley 
2003; Johansson and Ericson 2004; Helbig et al. 2005), 
and others have been shown to be relatively invariant 
at or below the species level (Wickliffe et al. 2003; 
Carroll and Bradley 2005). Alternatively, Fonseca et 
al, (in press) and Porter et al. (in review) found intron 
sequences useful for resolving relationships among 
closely related species of bats. 

The purpose of this study was to examine the 
phylogenetic utility of DNA sequence data from the 
2nd intron of the alcohol dehydrogenase gene 1 (Adh 1) 
obtained for 41 species in six genera of rodents. Acih 
genes, in mammalian systems, code for dimeric zinc- 
metallo enzymes that catalyze oxidization of alcohols 
to aldehydes or ketones (Szalai et al. 2002). Seven 
distinct classes are recognized in vertebrates (Duester 
etal. 1999; Dolney et al. 2001; Szalai et al. 2002), six 
of which are present in mammalian species (Duester 
et al. 1999; Szalai et al. 2002). In rodents, four Adh 
classes (I, 11, 111, and IV) have been reported (Crabb 
and Edenberg 1986; Zhang et al. 1987; Park and Plapp 
1991; Bradley et al. 1993, 1998; Duester et al. 1999; 
Szalai et al. 2002) and are controlled by six genes, 
Adh\,Adh2,Adh2,AdhA,Adh5a, andAdh5b (Duester et 
al. 1999; Dolney et al. 2001; Szalai et al. 2002). Adh 1 
is approximately 13 kilobase pairs in length (Crabb et 
al. 1989), contains nine exons (coding between 6 and 
87 amino acids) and eight introns (ranging from 91 bp 
to 3.6 kilobase pairs; Zhang et al. 1987). Hereafter, in 
referring to the 2nd intron, we follow the recommended 
nomenclature (Adh l -12^ of vertebrate alcohol dehydro¬ 
genases outlined in Duester et al. (1999). 
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Materials and Methods 


Specimens examined. —Specimens examined, 
GenBank accession numbers, and information associ¬ 
ated with museum vouchers are listed in Appendix I. 
Complete sequences for Adh\-\2 and cytochrome-6 
( Cyt-b ) were obtained for 41 individuals and 10 indi¬ 
viduals, respectively, and were deposited in GenBank. 
In addition, 31 Cyt-b sequences were retrieved from 
GenBank. Holochilus chacarius was designated as the 
outgroup taxon in phylogenetic analyses of Adh\-\2 
and Cyt-b data, as previous morphological and mo¬ 
lecular studies agree that Holochilus may be outside 
the remainder of taxa in this study (reviewed in Musser 
and Carleton 2005). Sequence data were evaluated 
and relationships were inferred among the 40 ingroup 
species (13 Neotoma , six Oryzomys, one Ototylomys , 
19 Pern my sens, and one Tylomys) representing the 
subfamilies Neotominae and Sigmodontinae (Musser 
and Carleton 2005). 

Adh 1-12 data. —Genomic DNA was isolated from 
approximately 0.1 g frozen muscle or liver tissues by the 
method of Smith and Patton (1999). Polymerase chain 
reaction (PCR; Saiki et al. 1988) primers for Adh 1-12 
(Table 1) were designed from conserved regions in 
exons 2 and 3 based on sequences from Homo (Ikuta 
et al. 1986), Geomys (Bradley et al. 1993, 1998), Mns 
(Zhang et al. 1987), Peromyscus (Zheng et al. 1993), 
and Rattus (Crabb et al. 1989). A nuclear fragment 
approximately 650 bp long, encompassing the entire 
Adh\-\2, was amplified in all six genera using two for¬ 
ward primers (2340-1 or EXON 1I-F) and two reverse 
primers (2340-11 or EXON III-R). PCR thermal profiles 


varied only slightly among genera examined: initial 
denaturation at 95°C for 2-10 min, followed by 25-30 
cycles of denaturation at 95°C for 30-60 sec, annealing 
ramped from 52-53°C down to 46-48°C and back up to 
52-53°C at a rate of 0.6°C/sec, extension at 73°C for 
1.5 min, and a final extension cycle of 73°C for 4 min. 
Ramping speed between all three phases of PCR was 
set at a rate of 1 °C/sec. Reaction concentrations (35 pi 
volume) included approximately 300 ng genomic DNA, 
0.07 mM dNTPs, 2.86 mM MgCl, 3.5 pi 10X buffer, 
0.286 pM primer, and 1.25-1.5 U enzyme (FailSafe 
PCR Enzyme Mix, Epicentre, Valencia, California). 
In some reactions, 1.5 U AmpliTaq Gold (PE Applied 
Biosystems, Foster City, California) was used with an 
initial denaturation time of 10 min. 

PCR products were purified using the QIAquick 
PCR purification kit (Qiagen, Valencia, California) and 
sequenced with Big-Dye version 3.1 chain terminators 
(Applied Biosystems, Inc., Foster City, California). 
Appropriate external primers and four internal prim¬ 
ers (350F, 350R, 41 OF, 41 OR; Table 1) were used to 
sequence each strand entirely. Thermal profile for 
cycle sequencing was modified from the manufacturer’s 
recommendation: 25-30 cycles of 95°C for 30 sec dena¬ 
turation, 50°C for 20 sec annealing, and 60°C for 3 min 
extension. Sequencing reactions were precipitated and 
concentrated with standard isopropanol methods, re¬ 
suspended in 15ml Hi-Di Formamide, and electropho- 
resed on an ABI31 00-Avant Genetic Analyzer (Applied 
Biosystems, Inc., Foster City, California). Sequencher 
3.0 software (Gene Codes, Ann Arbor, Michigan) was 


Table 1. Primer sequences used in the polymerase chain reaction (PCR) amplification and sequencing 
of the alcohol dehydrogenase locus. Sequences are given in a 5 ' to 3 ' orientation. 

Primer 

Sequence 

2340-1 (5’PCR primer) 

2340-11 (3’ PCR primer) 

EXON II-F (5’ PCR primer) 

EXON III-R (3’PCR primer) 

350F (internal sequencing primer) 

350R (internal sequencing primer) 

41 OF (internal sequencing primer) 

41 OR (internal sequencing primer) 

GTAATCAAGTGCAAAGCAGCTG 

TAACCACGTGGTCATCTGAGC’G 

GTAATCAAGTGCAAAGCRGCYYTRTGGGAG 

GACTTTATCACCTGGTTTYACWSAAGTCACCCC 

GTGCTAAACATCTTGATTCCRAAAG 

GCTTTTGGAATCAAGATGTTTAG 

CTATAGCACAGCACAGC 

TGCTGTGCTGTGCTATAG 
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used to assemble and proof resultant fragments. Base 
calling ambiguities on single strands were resolved 
by choosing the call on the cleanest strand or by us¬ 
ing appropriate IUB ambiguity code if both strands 
showed the same ambiguity. 

Because provisional statements of homology 
(i.e., sequence alignment) are of special concern for 
non-coding DNA sequences possibly containing in¬ 
sertion/deletions (Giribet and Wheeler 1999), and be¬ 
cause many parameters can affect multiple-sequence 
alignment and resulting phylogenetic inference (De- 
Salle et ah 1994; Hickson et al 2000; Lutzoni et al. 
2000; Wheeler 1995), two multiple-sequence align¬ 
ments o^Adh 1 -12 data were performed independently 
in Clustal X software (Thompson et al. 1997): one 
using a 30:4 gap cost-ratio, the other using a 5:4 gap 
cost-ratio (Van Den Bussche and Hoofer 2001; Van 
Den Bussche et al. 2002). Alignments subsequently 
were examined using MacClade software (version 
4.05; Maddison and Maddison 2002), ambiguously 
aligned sites were delimited following methods of 
Hoofer and Van Den Bussche (2003), and analyses 
were performed with and without those sites. 

Cyt-b data .—Sequences of the entire (1,143 
bp) mitochondrial Cyt-b gene were obtained from 
each of the 41 species to facilitate comparison to 
the Adh\-\2 data. PCR and sequencing primers, 
conditions, and thermal profiles followed Edwards 
and Bradley (2002) and Bradley et al. (2004a). 
Sequence alignments were performed manually 
and checked in MacClade software (version 4.05; 
Maddison and Maddison 2002) to ensure there were 
no insertions/deletions or stop codons in the protein¬ 
coding gene. 

Data analysis .—Analyses were performed 
in PAUP* 4.0b 10 software (Swofford 2002) or 
MrBayes 2.01 software (Huelsenbeck and Ronquist 
2001). Nucleotide positions were treated as unor¬ 
dered, discrete characters (A, C, G, and T), multiple 
states as polymorphisms, and gaps as missing. 
Nucleotide sequences from both Adh \-12 and Cyt- 
b were evaluated four ways: 1) base frequencies, 
number of transitions, number of transversions, and 
number of substitutions per 100 bp were estimated 
within and compared among the five ingroup genera; 
2) levels of phylogenetic signal were estimated using 


the ^ r 1 -statistic (Hillis and Huelsenbeck 1992) for 
100,000 randomly drawn trees; 3) genetic distances 
(uncorrected “p”) were obtained and compared using 
pairwise comparisons of taxa; and 4) maximum like¬ 
lihood, Bayesian likelihood, and parsimony analyses 
were performed and compared among taxa. 

Based on hierarchical likelihood ratio tests 
(hLRTs) in Modeltest software (Posada and Cran¬ 
dall 1998), the following models of nucleotide 
substitution and associated parameters best fit the 
data; Adh\-\2 —Hasegawa-Kishino-Yano (HKY) 
model with allowances for gamma distribution of rate 
variation (G), ti/tv = 2.308, JtA = 0.301, JtC = 0.191, 
jiG = 0.183, JtT = 0.325, a = 1.662; Cyt-b —general 
time reversible (GTR) with allowances for G and 
proportion of invariant sites (I), R A( = 1.310, R A( . = 
10.479, R at = 2.473, R C( . = 0.708, R CT = 32.137, jtA 
- 0.375, JtC = 0.343, jtG = 0.079, JtT = 0.203, a = 
0.687, p - 0.456, The HKY + G model best fit the 
Adh 1 -12 data with and without ambiguous characters, 
although specific model parameters differed slightly; 
Adh 1-12 values reported were calculated without 
ambiguous characters. 

Bayesian analyses were run at least two mil¬ 
lion generations with four Markov-chains, random 
starting trees for each chain, and trees sampled 
every 100th generation. For each data set, two 
independent analyses were run to assess whether 
chains converged on the same posterior probability 
distribution and whether likelihood values became 
stable (Huelsenbeck et al. 2002). Model parameters 
were treated as unknown variables (with uniform 
priors) to be estimated in each Bayesian analysis 
(Leache and Reeder 2002). Burn-in values (initial 
set of unstable generations to be ignored) were based 
on empirical evaluation of likelihoods converging 
on stable values. Clade reliabilities were assessed 
using posterior probabilities (values > 0.95 regarded 
as significant). 

Maximum likelihood analyses were performed 
with full heuristic searches, neighbor-joining start¬ 
ing trees, and tree-bisection-reconnection branch 
swapping. Parsimony analyses, with all characters 
and substitution types given equal probabilities (i.e., 
unweighted), were conducted with full heuristic 
searches with 10 random additions, starting trees 
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by simple addition, and tree-bisection-reconnection 
branch swapping. Clade reliabilities were assessed 
in parsimony analyses using bootstrapping methods 
with 250 iterations (Felsenstein 1985). Due to pro¬ 


hibitive computation time under maximum likelihood, 
bootstrapping was performed for 100 iterations for just 
the Adh 1-12 dataset (100 iterations). Values > 70 were 
regarded as strong support. 


Results 


Adh 1-12 data .—-Complete sequence of Adh 1- 
12 averaged 530 bp for the 41 rodents examined, 
ranging from 494 ( Holochilus) to 576 (Peromyscus 
cal if or mens). Alignment of sequences resulted in 614 
aligned sites (125 ambiguously aligned) with the 5:4 
gap cost-ratio, and 607 aligned sites (133 ambiguously 
aligned) with the 30:5 gap cost-ratio, corresponding 
to the insertion of more gaps with a lower cost ratio. 
As results from all subsequent analyses, phylogenetic 
analyses in particular, essentially were identical regard¬ 
less of alignment and with and without the ambiguously 
aligned characters excluded, only results based on the 
5:4 alignment with 125 ambiguous characters removed 
are reported and discussed. 

After removing 125 ambiguous characters, 489 
characters were available for analysis, of which 260 
were constant and 126 were parsimony informative. 
Overall nucleotide frequencies varied slightly among 
the three genera (Table 2), averaging 30.47% (A), 
17.51% (C), 17.94% (G), and 34.08% (T). The transi¬ 
tion to transversion ratio was approximately 2.28 to 1. 
The number of heterozygous sites ranged from zero 
(37 taxa) to three (N. Stephens /, N. mexicana , and O. 
perenensis ), with a mean heterozygosity of 0.34 per 
taxon. The g { statistic was skewed significantly left 


(-0.65; P < 0.01), indicating strong phylogenetic signal 
(Hillis and Huelsenbeck 1992). 

Cyt-b data. —Complete sequences of Cyt-h (1,143 
bp) were included for the 41 rodents examined herein. 
Sequence alignment was unequivocal and contained no 
stop codons. Of the 1,143 characters, 589 were constant 
and 475 parsimony informative. Nucleotide variation 
was distributed across codon positions as expected for 
protein-coding genes (1st position, 111; 2nd position, 
37; 3rd position, 327). Overall nucleotide frequen¬ 
cies varied slightly among the three genera (Table 2), 
averaging 31.67% (A), 28.15% (C), 12.45% (G), and 
27.74% (T), and transition to transversion ratio was 
approximately 2.68 to 1. The g ] statistic was skewed 
significantly left (-0.49; P < 0.01), indicating strong 
phylogenetic signal (Hillis and Huelsenbeck 1992). 

Phylogenetic analyses. —For Adh 1 -12 and Cyt-b , 
Bayesian likelihoods reached stability before 100,000 
generations (i.e., bum-in = 1,000), thinning the data to 
19,000 sample points. Topology and posterior prob¬ 
abilities for nodes and model parameters for all sets of 
runs agreed. Maximum likelihood analysis resulted in 
a single best tree for both Adh 1-12 (Lnl = -2,718.63) 
and Cyt-b (Lnl = -14,379.92) data sets (Figs. 1 and 2). 


Table 2. Nucleotide base composition (A, C, G, and T), number of transition substitutions per 100 bp (Wi/100 bp), 
number of transversion substitutions per 100 bp (#Tv/l 00 bp), and the uncorrected p distance (p Dist) for the 3 genera 
possessing multiple taxa. All values were averaged across taxa within each of the 3 genera for nitron 2 of the alcohol 
dehydrogenase locus (Adh 1-12) and the mitochondrial cytochrome-b gene (Cyt-b). 


Taxon/Sequence 

% A 

%C 

% G 

% T 

# Ti/100 bp 

# Tv/100 bp 

p Dist 

Peromyscus Aclh\-\2 

30.16 

17.48 

18.58 

33.77 

2.31 

1.09 

3.41% 

Peromyscus Cyt-b 

31.72 

27.45 

12.62 

28.21 

8.38 

2.79 

12.29% 

Neotoma Adh 1 -12 

30.35 

17.73 

17.59 

34.33 

1.67 

0.50 

2.22% 

Neotoma Cyt-b 

32.24 

28.93 

12.43 

26.40 

8.73 

2.90 

11.61% 

Oryzomys Adh 1 -12 

31.34 

17.33 

16.97 

33.37 

4.48 

2.10 

6.70% 

Oryzomys Cyt-b 

31.16 

28.43 

11.98 

28.43 

7.73 

5.75 

14.00% 
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Parsimony analysis resulted in 62 most-parsimonious 
trees (length = 363, Cl = 0.78, Ri = 0.89) and a single 
most-parsimonious tree (length = 3,321, Cl = 0.28, Rl 
= 0.48) for Adh 1-12 and Cyt-b , respectively. 

There were some topological differences within 
and between datasets and between the three optimality 
criteria, although only one difference between datasets 
represented a statistically supported conflict (i.e., > 70% 
bootstrap value, >0.95 Bayesian posterior probability). 


This conflict involved the sister species relationships 
of N. albigula , N. leu codon, and N. micropus. Adh 1-12 
supported a N. cilbigula/N. leucodon sister relationship, 
whereas Cyt-b supported a N. leucodoniN . micropus 
sister relationship (Figs. 2 and 3). Overall, statistically 
supported topologies obtained from all optimality cri¬ 
teria agreed within and between datasets, supporting 
monophyly of Neotoma and Peromyscus and diphyly 
of Oryzomys (Fig. 3). 


Discussion 


Results of this study indicated that Adh 1 -12 DN A 
sequences are useful for addressing phylogenetic re¬ 
lationships within and among the Sigmodontinae and 
Neotominae. Based on the g statistic and consistency 
and retention indices, the Adh 1 -12 dataset had a greater 
overall phylogenetic signal and less homoplasy than 
the Cyt-b dataset. Adh 1-12 sequences also provided 
statistically supported resolution to slightly more re¬ 
lationships than Cyt-b (Fig. 3). For example, Adh\-\2 
strongly supported all relationships examined at and 
above the genus level, including monophyly of Neo¬ 
toma and Peromyscus and a sister relationship between 
Ototylomys and Tylomys , relationships that have been 
documented repeatedly by analysis of morphological 
and DN A sequence data (Carleton 1980; Bradley et al. 
2004b; Edwards and Bradley 2002; Reeder and Bradley 
2004, in press a and b). Adh 1 -12 supported a diphyletic 
origin for Otyzomys , a relationship also supported by 
Cyt-b and recent studies of nuclear markers (Myers and 
Tucker 1995; Smith and Patton 1993; Weksler 2003), 

It is particularly noteworthy that, compared with 
Cyt-b , the Adh 1-12 dataset provided supported resolu¬ 
tion to more species-level relationships within Neoto¬ 
ma , Peromyscus , and Oryzomys (Fig. 3), despite fewer 
overall characters (489 Adh 1 -12, 1,143 Cyt-b) and fewer 
parsimony informative characters (126 Adh 1-12, 464 
Cyt-b). Furthermore, all but one of the species-level 
relationships were congruent with results from Cyt-b. 
The exception involved the sister species relationships 
of A. albigula , N. leucodon , and N. micropus , and per¬ 
haps represented a situation of misleading phylogenetic 
inference. Whereas Adh 1-12 analysis provided strong 
support for a sister relationship between N. albigula 


and N. leucodon (parsimony, 98%; maximum likeli¬ 
hood, 94%; Bayesian, 100%), Cyt-b analysis strongly 
supported a sister relationship between A. leucodon and 
N. micropus (parsimony, 86%; Bayesian, 100%). The 
difference, however, is that Adh 1-12 support was based 
on three characters, whereas Cyt-b support was based 
on 22 characters. We view these particular Adh 1-12 
results with caution, and suggest additional study of 
nuclear markers will be necessary to help determine 
whether Adh 1-12 results accurately reflect a conflict 
between the mitochondrial and nuclear genomes or 
represents the case of misleading inference of relation¬ 
ship due to inadequate variability. 

Overall, the Adh 1-12 results are congruent with 
results from Cyt-b and other types of data, and gener¬ 
ally are accompanied by high statistical support. Fur¬ 
thermore, multiple sequence alignment, which can be 
problematic with non-coding DNA sequences (Giribet 
and Wheeler 1999; Hoofer and Van Den Bussche 2003; 
Van Den Bussche et al. 2002), was not of particular 
concern in this study. Although we identified 125 
ambiguously aligned characters in the 5:4 alignment 
and 133 in the 30:5 alignment, results from analysis 
of both alignments with and without the ambiguous 
characters were identical in topology and statistical 
support. Thus, most of the phylogenetic signal in the 
Adh\-\2 dataset was associated with polymorphisms 
(nucleotide substitutions) rather than insertions/dele¬ 
tions. Yet, some insertion/deletions were informative. 
For example, a 5 bp (nucleotide positions 335-339) and 
a 4 bp (nucleotide positions 487-490) insertion/deletion 
were present in species of Oryzomys and Holochilus 
but absent in Neotoma and Peromyscus. Similarly, a 17 
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— 0.01 substitutions/site 


P, attwateri 
— P. difficilis 
L P. grafrvs 
|—P. melanophrys 
1 — P. perfulvus 
*— P. megalops 
P. beatae 
P. boylii 

- P. hylocetes 

- P. levipes 

- P. schmidlyi 
L P. spicilegus 

P. californicus 

-P. eremicus 

— P. leucopus 
-P. maniculatus 


— P crinitus 
— P. mexicanus 
— P. pectoralis 
r N. leucodon 
'— N. albigula 
— N. mexicana 
— N. goldmani 
j— N. floridana 
I_' N. magister 


— N. micropus 
r— N. picta 
‘ N. isthmica 
N. stephensi 

N. fuscipes 
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N. cinerea 
*— N. lepida 


Ot. phyllotis 
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■ O. albigularis 
-O. perenensis 

— O. melanotis 

— O. alfaroi 


■ O. couesi 
O. palustris 
■ Holochilus 


Figure 1, Maximum likelihood phylogram (Lnl -2,718.627) from analysis of 5:4 gap cost-ratio alignment 
of complete Adh 1-12 sequences (489 base pairs; 614 aligned sites minus 125 ambiguously aligned sites) 
using best-fit model (HKY + G; ti/tv = 2.308, jtA = 0.301, jtC = 0.191, jiG = 0.183, jtT = 0.325, a = 
1.662). Holochilus was the designated outgroup. /V. = Neotoma. O. = Oryzomys. Ot. = Ototylomys. 
P. = Peromyscus, T. = Tylomys. 
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. — N. magister 
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N. leucodon 
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r~ N. picta 
■ N. isthmica 
— N. mexicana 
- N. cinerea 


— N. stephensi 

- N. fuscipes 

N. lepida 


T. 


j- P. leucopus 

'— P. maniculatus 
r“ P. melanophrys 
P. perfulvus 
— P. mexicanus 




P. megalops 
P. attwateri 
P. difficilis 

-P. pectoralis 

r P beatae 
r L P. levipes 

- L P. schmidlyi 
. — P. boylii 
_p P. hylocetes 

— P. spicilegus 

— P. gratus 

~ P californicus 

■ P eremicus 

■ P crinitus 


- Ot. phyllotis 
nudicaudus 


— O. alfaroi 

— O, melanotis 
O. perenensis 




— O. albigularis 
O. couesi 
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Figure 2. Maximum likelihood phylogram (Lnl -14,379.922) from analysis of complete Cyt-b 
sequences (1,143 base pairs) using best-fit model (GTR + G + I; R AC = 1.310, R APi = 10.479, R A1 
- 2.473, R cg = 0.708, R CT - 32.137, jrA= 0.375, jrC - 0.343, jtG - 0.079, JtT = 0.203, a = 0.687, 
p = 0.456). Holochilus was the designated outgroup. N. = Ncotomu, O. ~ Oryzomys. Ot . = 
Ototylomys. P. = Peromyscus . T. — Tylomys. 
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Figure 3. Cladograms from phylogenetic analyses of complete Adh 1-12 (left) and Cyt-b (right) DNA sequences, 
Holochilus was the designated outgroup. In order, numbers above and below nodes are bootstrap proportions (250 
iterations) from parsimony analysis, bootstrap proportions (100 iterations) from maximum likelihood analysis, and 
posterior probability proportions from Bayesian analysis. Maximum likelihood bootstrapping was feasible only for 
Adh\-\2 data. Only nodes with >70% bootstrap support or >0.95 posterior probabilities, or both, are shown. Vertical 
dotted line indicates sister relationship between TV. micropus and clade of TV. floridana + N. magister. N. = Neotoma. 
O. - Oiyzomys. Ot. = Ototylomys. P. - Peromyscus. T, = Tylomys. 
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bp insertion/deletion (nucleotide positions 502-518) 
was present in Neotoma and Peromyscns but was ab¬ 
sent in Oryzomys and Holochilus . These results are 
encouraging toward the goal of recovering reliable 
phylogenetic relationships at low and intermediate 


taxonomic levels from a nuclear intron. Further 
study of this intron, along with other nuclear and 
mitochondrial markers, should aid our understand¬ 
ing of the phyletic limitations oi Adh\-Yl as well as 
organismal genealogy. 
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Appendix I 

Specimens examined. — Collection localities, museum acronyms, and GenBank accession numbers are 
provided for each specimen examined in this study. Specimens are from the United States unless otherwise noted. 
Abbreviations for museum acronyms (in parentheses and to the left of the semicolon) follow Hafner et al. (1997): 
Abilene Christian University Natural History Collection (ACUNHC), Angelo State University Museum Natural 
History Collections (ASNHC), Monte L. Bean Life Science Museum (BYU), Texas Cooperative Wildlife Col¬ 
lection (TCWC), Museum of Texas Tech University (TTU), The Museum of Southwestern Biology (MSB), and 
UniversidadNacional Autonoma de Mexico (UNAM). If museum catalogue numbers were unavailable, specimens 
were referenced with the corresponding TK number ( special number of the Museum of Texas Tech University). 
GenBank accession numbers (AF, AY, U, and DQ) for Adh\-\2 and Cyt-b are provided in parentheses to the right 
of the semicolon and separated by a comma, respectively. 

Holochilus chacarins. —PARAGUAY: Departmento Pte. Hayes; Estancia Loma Pora, 23°33.15’S, 57°34.3’W (TTU104423; 
DQ227456, DQ227455). 

Neotoma albigula. —New Mexico; Yuma Co., 3.7 km S, 5.6 km W Somerton, UTM ! 1 708569E 3608362N (TTU78451; AY817648, 
AF376477). 

Neotoma cinerea. —Utah; San Juan Co., Owaehamo Bridge (MSB 121427; AY817635, AF 186799). 

Neotoma fioridana. —South Carolina; Richland Co.. Congaree Swamp National Monument, 33°49’N, 80°50’W (MSB74955; 
AY817637, AF294335). 

Neotoma juscipes. —California; Riverside Co., Rancho Capistrano (Ortega Mountains) (TTU81391; AY817632, AF376479). 

Neotoma golduumi, —MEXICO: Nuevo Leon; 1 km S Providencia (TTU45227; AY817656, AF 186829). 

Neotoma isthmica.— MEXICO: Oaxaca; Las Minas, UTM 15 191165E 1824954N (TTU82665; AY817630, AF329079). 

Neotoma lepida. —California; Orange Co., Irvine Lake, 1.3 km E Fremont Canyon on Lake View access road (TTU79131; AY817633, 
AF307835). 

Neotoma leucodon. —Texas; Kerr Co., Kerr Wildlife Management Area, UTM 14 452336E 3330772N (TTU71198; AY817643, 
AF 186806). 

Neotoma magister .— Virginia; Madison Co., Shenendoah National Park, White Oak Canyon, 38°34 , 36”N, 78°22’30”W (MSB74952; 
AY817641, AF294336). 

Neotoma mexicana .— Texas; Jeff Davis Co., Mount Livermore Preserve, UTM 13 579953E 3389871N (TTU101643; AY817645, 
AF294346). 

Neotoma micropus .— New Mexico; Roosevelt Co., 26.4 km S, 4.8 km E Taiban (TTU catalogue number unavailable, TK31643; 
AY817652, AF 186822). 

Neotomapicta .— MEXICO: Guerrero; 6.4 km SSW Filo de Caballo (UNAM catalogue number unavailable, TK93390; AY817629, 
AF305569). 

Neotoma stephensi .— Arizona; Navaho Co., 4.8 km S Woodruff, UTM 12 588361E 3844338N (TTU78505; AY817642, 
AF308867). 

Oryzomys albigitlaris .— ECUADOR: Pichincha; 10 km NW Quito, Tandayapa Valley, 0°00T3”N, 78°40’70’AV (ACUNHC917; 
DQ207945, DQ224407). 

Oiyzomys alfaroi .— NICARAGUA: Selava Negra; Altajo Trail (TTU 101644; DQ207950, DQ224410). 
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Appendix 1 (cont.) 


Oryzomys cousei. —HONDURAS: Olancho; 4 km E Catacamas (Escuela de Sembrador), UTM 16 624523E 163751 IN (TTU84697; 
DQ207948, DQ185383). 

Oryzomys melanotis. —HONDURAS: Atlantida; Lancetilla Botanical Gardens, UTM 16 451012E 1740282N (TTU84374; DQ207947, 
DQ224409). 

Oryzomyspalnstris. —Texas; Galveston Co., Texas City, Virginia Point (TTU82920; DQ207949, DQ 185382). 

Oryzomys perenensis — PERU: Loreto; Maynas, 25 km S Iquitos (Estacio Biologia Allpahuayo) (TTU98606; DQ207946, 
DQ224408). 

Ototylomysphyllotis. — HONDURAS: Atlantida; Lancetilla Botanical Gardens, UTM 16 451012E 1740282N (TTU84371; AY817624, 
no Cyt-b data); MEXCIO: Quintana Roo, 1 km N Noh-bec (ASNHC7254; no Adh\-\2 data. AY009788). 

Peromyscus attwateri —Oklahoma; McIntosh Co., 5 km E Dustin (TTU55688; AY817626, AF155384). 

Peromyscus beatae. —MEXICO: Chiapas, Yalentay, UTM 15 52417 IE 1852486N (UNAM catalogue number unavailable, TK93279; 
AY994223, no Cyt-b data); Veracruz; Xometla (TCWC48060; no Adh\-\2 data, AF131921). 

Peromyscus boylii. —California; San Diego Co., Heise County Park (TTU83102; AY994225, no Cyt-b data); MEXICO: Aguascali- 
entes; Rincon de Romos (TCWC48438; no Adh\-\2 data, AF131924). 

Peromyscus californicus. —California; Los Angeles Co., Chatsworth Resevoir Park (TTU83292; AY994211, AF155393). 

Peromyscus crinitus. —Utah; Emery Co., Cottonwood Canyon, 39°16’5 1.8”N, 111°10’31.9”W (BYU18639; AY994213, 
AY376413). 

Peromyscus difficilis. —-MEXICO: Tlaxcala; 2 km NE Tepetitla (TTU82690; AY994219, AY376416). 

Peromyscus eremicus, —California: Los Angeles Co., Kanan Dome Road (TTU81850; AY994212, no Cyt-b data); California; Los 
Angeles Co., Calabasas Creekside Park (TTU83249; no Adh\-\2, AY322503). 

Peromyscus gratus.— MEXICO: Michoacan; 4 km E Costzeo (UNAM catalogue number unavailable, TK46354; AY994218, 
AY376421). 

Peromyscus hylocetes. —MEXICO: Michoacan; Estacion Cerra Burror, Microodas; 3,270 m (UNAM catalogue number unavailable, 
TK 45309; AY994235, DQ000481). 

Peromyscus leucopus .—Texas; Presidio Co., Las Palomas Wildlife Management Area, UTM 13 52948 IE 332I292N (TTU75694; 
AY994240, no Cyt-b data); Dickens Co., 0.9 km E Afton (TTU catalogue number unavailable, TK47506; no Adh\- 12 data, 
AF131926). 

Peromyscus levipes ,—MEXICO: Michoacan, Las Minas, 3 km SW Tuxpan (UNAM catalogue number unavailable, TK.47819; 
AY994224, DQ000477). 

Peromyscus maniculatus. —Arkansas; Mississippi Co., Dillahunty Pecan Orchard (TTU97830; AY994242, no Cyt-b data); Alaska, 
Alexander Archipelago, Bushy Island, Petersburg Quad, 56°15’45”N, 132°58’52”W (UAM50770; noAdh\-\2 data, AF119261). 

Peromyscus megalops. —MEXICO: Guerrero; 6.4 km SSW Filo de Caballo (TTU82712; AY994217, DQ000475). 

Peromyscus melanophrys. —MEXICO: Durango; 2.2 km S, 2.5 km E Vicente Guerrero (TTU75509; AY9942I6, AY322510). 

Peromyscus mexicanus. —MEXICO; Chiapas; 14.4 km N Ocozocoaulta (TTU82759; AY994236, AY376425). 





Appendix l (cont.) 


Peromyscus pectoralis .—MEXICO: Jalisco; 30 km W Huejuquilla del Alto (UNAM catalogue number unavailable, TK48645; 
AY994221, no Cyt-b data; TTU75575; no Adh\M data, DQ000476). 

Peromyscus perfulvus .—MEXICO: Michoacan; Tunel de Riego, 2 km E Cerro Colorado, 1290 M, 19°19’220”N, 100°28’308”W 
(UNAM catalogue number unavailable, TK47926; AY994215, DQ000474). 

Peromyscus schmidiyi —MEXICO: Durango; 6.1 km W Coyotes, Hacienda Coyotes, UTM 13 465908E 2634281N (TTU81617; 
AY994228, AY322524). 

Peromyscus spicilegus .—MEXICO: Michoacan; KM 81 carr., Ario de Rosales-La Huacana, 1602 m, 19 o 10’59”N, 101°43’42 ,, W 
(UNAM catalogue number unavailable, TK 47888; AY994232, DQ000480). 

Tylomys midicaudatus .—GUATEMALA: lzabal, Cerro San Gil (TTU62082; AYS 17625, AF307839). 
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